Ab initio many-body calculation of the "^66(^,7)^6 radiative capture 
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We apply the ab initio no-core shell model/resonating group method (NCSM/RGM) approach 
to calculate the cross section of the ''Be(p,7)*B radiative capture. This reaction is important for 
understanding the solar neutrino flux. Starting from a selected similarity-transformed chiral nucleon- 
nucleon interaction that accurately describes two-nucleon data, we performed parameter-free many- 
body calculations that simultaneously predict both the normalization and the shape of the S-factor. 
We study the dependence on the number of ^Be eigenstates included in the coupled-channel equa- 
tions and on the size of the harmonic oscillator basis used for the expansion of the eigenstates and 
of the localized parts of the integration kernels. Our S-factor result at zero energy is on the lower 
side of, but consistent with, the latest evaluation. 



The core temperature of the Sun can be determined 
with high accuracy through measurements of the *B neu- 
trino flux, currently known with a 9% precision [l|. An 
important input in modeling this flux is the ^Be(p,7)®B 
reaction Q that constitutes the final step of the nucle- 
osynthetic chain leading to ®B. At solar energies this 
reaction proceeds by external, predominantly nonreso- 
nant El, S- and D-wave capture into the weakly-bound 
ground state (g.s.) of ®B. Experimental determinations 
of the ''Be(p,7)*B capture include direct measurements 
with proton beams on ^Be targets as well as indi- 

rect measurements through the breakup of a ^B projec- 
tile into ''Be and proton in the Coulomb field of a heavy 
target Theoretical calculations needed to extrapo- 

late the measured S-factor to the astrophysically relevant 
Gamow energy were performed with several methods: the 
R- matrix parametrization Q , the potential model [lol - 
[l3 |. microscopic cluster models [l3- 15| and, recently, also 
using the ab initio no-core shell model wave functions for 
the ^B bound state The most recent evaluation 

of the ^Be{p,7)®B S-factor at zero energy, £'17(0), has a 
~10% error dominated by the uncertainty in theory 

In this Letter, we present the first parameter-free ab 
initio many-body calculations of the ''Be(p,7)®B capture 
starting from a nucleon-nucleon (NN) interaction that 
describes two-nucleon properties with high accuracy. We 
apply a recently developed technique that combines ab 
initio no-core shell model (NCSM) [l3| and resonating- 
group method (RGM) into a new many- body ap- 

proach [iO-il (ab initio NCSM/RGM) capable of treat- 
ing bound and scattering states of light nuclei in a unified 
formalism. We use, in particular, the orthonormalized 
NCSM/RGM many-body wave functions given by 



l*-^"^)=5] J drr^ J dr'r'^ A,\^ 



(1) 



mass separation r^^a.a, and binary-cluster channel states 

(.rX) - rA-a,a) 



X Yi {fA^a,a) 



rVA-a.a 



(2) 



The wave functions xi"^^ i''^) of the relative inter-cluster 
motion satisfy the integro-differential coupled-channel 
equations 



(3) 



with bound- or scattering-state boundary conditions. 
The Hamiltonian and norm kernels, 



Hf/(r',r) = ($fJ 



vr 



(4) 
(5) 



with the inter-cluster antisymmetrizer Ay, the center-of- 



contain all the nuclear structure and antisymmetrization 
properties of the problem. Further relevant details of the 
NCSM/RGM formalism are given in Ref. 20]. 

In the present case A is equal to 8, and the projectile 
is a proton \a—\ in Eq. ([5])]. The input into Eq. ([3]) are: 
(i) the chiral N'^LO NN potential [23|, which we soften 
by a similarity renormalization group (SRG) transforma- 
tion [H, characterized by an evolution parameter A; 
(ji) the eigenstates of the target, i.e. ^Be, calculated 
within the NCSM. In Fig. [1] we show the energy depen- 
dence of the ^Be g.s. on the harmonic-oscillator (HO) 
frequency (a) for the HO basis sizes A'^max = 4 to 12, 
with the VlWt results obtained using the importance- 
truncation scheme (26j . The frequency dependence is 
quite flat and, with the selected A^A^ potential, we reach 
converge for the g.s. at A^max ~ 12. The g.s. energy min- 
imum is found at h£l—VS> MeV and we choose this fre- 
quency for all subsequent calculations (including eigen- 
states and integration kernels). The convergence of the 
absolute energies of the lowest five ^Be states is presented 
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FIG. 1. (Color online) Calculated NCSM '^Be eigenenergies 
using the SRG-N^LO A^iV potential with A = 1.86 fm'^ . 
Panel (a) shows the dependence of the g.s. energy on the HO 
frequency for A'max ~ 4—12 (with A^max ~ 12 results obtained 
within the importance-truncated basis). Absolute energies of 
the lowest 5 eigenstates for A^'max = 4 — 10 and ftil = 18 MeV 
are compared to experimental values in panel (6). 



TABLE I. ^Be g.s. energy (in MeV), charge radius (in fm), 
g.s. quadrupole (in e fm^) and magnetic (in fijq) moments and 
Ml transition (in /x^) obtained within the NCSM using the 
SRG-N^LO NN potential with A = 1.86 fm~^ Experimental 
values are from Refs. [23, 1231 • 





Q 




B(Ml;l-^f" 


NCSM -38.46 2.46(2) 


-5.39(10) 


-1.15 


3.14 


Expt. -37.60 2.647(17) 




-1.3995(5) 


3.71(48) 



in panel (6) of Fig. [T] Compared to the experimental 
values, we observe a slight overbinding of the g.s. and 
an overestimation of the 7/2~ and 5/2^ state excitation 
energies, but, overall, the agreement is reasonable. In 
Table m we compare some of our (IT-)NCSM ^Be results 
based on calculations up to iVmax— 14 to experimental 
values. 

Using the five lowest A^max=10 eigenstates of ^Be, we 
first solve Eq. ^ with bound-state boundary conditions 




— 3/2 s=l 1=1 

3/2" s=2 1=1 

-- l/2"s=ll=l 

7/2" s=3 1=1 

- 5/2'| s=2 1=1 

5/2' s=2 1=1 
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FIG. 2. (Color online) Dominant P-wave components of the 
2+ *B g.s. wave function for iV^ax = 10 and hO. = 18 MeV, 
using the SRG-N^LO iVA^ potential with A = 1.86 fm"^ The 
NCSM/RGM calculation includes ^Be g.s. and 1/2", 7/2", 
5/2^ and 5/2^ excited states. The calculated s.e. is 136 keV. 



to find the g.s. of ^B. We note that the same iVmax (^max 
-|-1 for the positive parity states) value is used to expand 
the localized parts of the integrations kernels ^ and 
(P. The chosen SRG-N^LO NN potential with A=1.86 
fm~^ leads to a single bound state, 2+, with separation 
energy (s.e.) 136 keV, quite close to the experimental 
137 keV. For the calculation of the low-energy behavior 
of the 5*17 S-factor, a correct s.e. is very important. The 
fact that the experimental s.e. of ®B can be found us- 
ing the SRG potential with a A from a "natural" range, 
i.e. «1.8— 2.1 fm"^, is reassuring. In Fig. we plot the 
most significant components of the radial wave functions 
x{r) for the 2+ g.s. of *B. The dominant component is 
clearly the channel-spin s—2 P-wave of the ''Be(g.s.)-p 
that extends to a distance far beyond the plotted range. 
Remarkably, we notice a substantial contribution from 
the ^Be(5/22")-p P -wave. Clearly, for a realistic descrip- 
tion of the *B g.s., this state must be taken into account. 
The influence of still higher ^Be resonances on the S- 
factor results will be discussed at the end of this Letter. 

Next, we solve the same NCSM/RGM equations ^ 
with scattering-state boundary conditions for a chosen 
range of energies and obtain scattering wave functions 
and the scattering matrix. The resulting phase shifts and 
cross sections are displayed in Fig. [31 All energies are in 
the center of mass (cm.). We find several P-wave reso- 
nances in the considered energy range. The first 1"*" res- 
onance, manifested in both the elastic and inelastic cross 
sections, corresponds to the experimental ®B 1+ state at 

MeV (0.63 MeV above the p-^Be threshold) 
The 3''" resonance, responsible for the peak in the elas- 
tic cross section, corresponds to the experimental ^B 3"*" 
state at Ex=2.32 MeV. However, we also find a low-lying 
0+ and additional 1+ and 2+ resonances that can be dis- 
tinguished in the inelastic cross section. In particular, the 
s—1 P-wave 2+ resonance is clearly visible. There is also 
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FIG. 3. (Color online) Selected P-wave diagonal phase shifts 
of p-^Be elastic scattering (a), inelastic ^Be(p,p')'^Be(l/2~) 
cross section (6) and elastic '^Be(p,p)^Be differential cross sec- 
tion at Oc.m. = 148° (c). Calculations as described in Fig. [2] 



an s=2 P-wave 2+ resonance with some impact on the 
elastic cross section. These resonances are not included 
in the current A=8 evaluation [2^. We note, however, 
that the authors of the recent Ref. [s^ do claim obser- 
vation of low-lying 0^ and 2'^ resonances based on an 
R-matrix analysis of their p-''Be scattering experiment. 
Their suggested 0+ resonance at 1.9 MeV is quite close 
to the calculated O"*" energy of the present work. 

With the resulting bound- and scattering-state wave 
functions that are properly orthonormalized and anti- 
symmetrized ([T|), we calculate the ''Be(p,7)*B radiative 
capture using a one-body El transition operator. The 
resulting Sn factor is compared to several experimental 
data sets in panel (a) of Fig. SI In the data, one can see 
also the contribution from the 1+ resonance due to the 
Ml capture. Our calculated S-factor is somewhat lower 
than the recent Junghans data [5| but the shape repro- 
duces closely the trend of the GSI data Q and is quite 
similar to that obtained within the microscopic cluster 
model [3] used in the most recent 6*17 evaluation Q. 
The contributions from the initial 1~, 2~ and 3~ partial 
waves are shown in panel (b) of Fig. 01 Our calculated 
S'i7(0)wl9.4 eV b is on the lower side, but consistent with 




FIG. 4. (Color online) Calculated ■^Be(p,7)*B S-factor as 
function of the energy in the cm. compared to data (a) . Only 
El transitions were considered. Initial-state partial wave con- 
tributions are shown in panel (fe) . Calculation as in Fig. [2] 



the latest evaluation 20.8±0.7(expt)±1.4(theory) eV b. 

We studied the convergence of the ^Be NCSM calcu- 
lations in Fig. [T] To verify the behavior of our S-factor 
with respect to HO basis size and number of included 
^Be eigenstates, we performed additional calculations as 
summarized in Fig. [S] To study the dependence on the 
HO basis size, shown up to 7Vmax=12 in panel (6), we 
use the importance truncation scheme and, due to com- 
putational limitations, we include only the three lowest 
eigenstates of ''Be. The iVmax=10 and 12 S-factors are 
very close. In panel (a), we present results with up to 
8 ''Be eigenstates obtained in a iVmax=8 basis. Calcula- 
tions with more than 5 eigenstates are presently out of 
reach for larger A'max values. We can see a significant 
impact of the 5/2~ states (with only three ^Be states, 
*B is unbound in this case). Among the others only the 
8th state, 7/2^, contributes somewhat to the s.e. and 
fiattens the S-factor at higher energies. We note that 
we selected different SRG-N^LO NN potentials with the 
aim to match closely the experimental s.e. in each of the 
largest calculation. From these results we conclude that 
the use of the A'max^lO space is justified and a limitation 
to the five lowest ^Be eigenstates is quite reasonable (or 
that the A',„ax=8 space is insufficient and a limitation to 
just 3 states is unrealistic). Also, based on these results 
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FIG. 5. (Color online) Convergence of the ^Be(p,7)*B S-factor 
with the number of ^Be eigenstates (a) and the size of the HO 
basis used to expand the ^Be eigenstates and localized parts 
of the integration kernels (6) . The number of eigenstates and 
the calculated separation energy in each case is shown in the 
legend. HO frequencies of Ml = 19 MeV (a) and 17 MeV (b) 
corresponding to the respective minima of ■'Be g.s. were used. 



we estimate the uncertainty of our calculated 517(0) to 
be ±0.7 eV b. 

In conclusion, we performed parameter-free ab initio 
many-body calculations of the ^Be(p,7)^B radiative cap- 
ture that predict simultaneously both the normalization 
and the shape of the S-factor. Our S-factor result at 
zero energy, 5'i7(0)=19.4(7) eV b, is on the lo"wer side of, 
but consistent "with, the latest evaluation, and its shape 
follo"ws closely the data from Ref. Our calculations 
can be further improved by including effects of additional 
higher-lying ^Be resonances. This can be best done by 
coupling the NCSM/RGM binary-cluster basis -with the 
NCSM calculations for ^B as outlined in Ref. [3l|. The 
inclusion of three-nucleon interactions, both chiral and 
SRG-induced [s^l, is also desirable. Efforts in these di- 
rections are under "way. 
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